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ABSTRACT 

We present a new method based on phase analysis for the Galaxy and foreground component separation 
from the cosmic microwave background (CMB) signal. This method is based on a prevailing assumption 
that the phases of the underlying CMB signal should have no or little correlation with those of the 
foregrounds. This method takes into consideration all the phases of each multipole mode (I < 50, 
—£< m < £) from the whole sky without galactic cut, masks or any dissection of the whole sky into 
disjoint regions. We use cross correlation of the phases to illustrate that significant correlations of the 
foregrounds manifest themselves in the phases of the WMAP 5 frequency bands, which are used for 
separation of the CMB from the signals. Our final phase-cleaned CMB map has the angular power 
spectrum in agreement with both the WMAP result and that from Tegmark, de Oliveira-Costa and 
Hamilton (TOH), the phases of our derived CMB signal, however, are slightly different from those of the 
WMAP Internal Linear Combination map and the TOH map. 



1. INTRODUCTION 

The release of the first year results from the Wilkin- 
son Microwave Anisotropy Probe (WMAP ) opens a new 
epoch of the full-sky CMB investigation, particularly for 
the analysis of different foreground components (Ben- 

■ nett et al. 2003a, b,c; Hinshaw et al. 2003a, b). The analy- 
ses from different groups, apart from the WMAP results, 
have raised two main issues: firstly, why is the power 
of the quadrupole component of the CMB signal con- 
siderably suppressed? (Tegmark, de Oliveira-Costa & 
Hamilton 2003; de Oliveira-Costa, Tegmark, Zaldarriaga 
& Hamilton 2003; Efstathiou 2003) Secondly, is the CMB 
signal Gaussian? (Komatsu et al. 2003; Chiang, Nasel- 

■ sky, Verkhodanov & Way 2003; Park 2003; Eriksen et al. 
2003) The key point to these questions may lie in the area 
of component separation, related with the removal of the 
foregrounds from the all frequency maps. 

The WMAP group produced an Internal Linear Combi- 
nation (ILC) map (Bennett et al. 2003c) using a weighted 
combination of the 5 bands outside the Galactic plane. 
Tegmark, de Oliveira-Costa & Hamilton (2003) (hereafter 
TOH) perform an independent foreground analysis from 
the WMAP data using weighting coefficients dependent 
not only on the angular scales (£) but also on the Galactic 
latitudes. 

In this paper we propose an alternative method with 
the assumption that the phases of the CMB signal derived 
from foreground cleaning should correlate as little as pos- 
sible with those of the foregrounds. According to the sim- 
plest inflation theory (see for review Komatsu et al. 2003), 
pure CMB signals constitute a Gaussian random field, 
which have random and uniformly distributed phases. 



Because of the pronounced non-Gaussianity of the fore- 
grounds, it therefore becomes straightforward to regard 
the non-Gaussianity, if detected, from the WMAP Internal 
Linear Combination (ILC) map, the TOH maps or any 
other derived maps from the WMAP data as contami- 
nations of the foregrounds at various levels. As such, 
we should have significant cross-correlations between the 
phases of the foregrounds and the CMB signal. 

We perform the phase analysis of the signals for all 
the K-W bands of the WMAP in order to remove the 
Galaxy and extragalactic foregrounds. Below we will call 
our derived map the Phase-Cleaned Map (PCM), and our 
method the PCM method. Different from other meth- 
ods for component separation, our method is based sig- 
nificantly on the properties of phases, using each map 
from the WMAP K-W bands as an image without any 
assumptions about the statistical nature of the CMB sig- 
nal. Moreover, we use the whole sky for analyses and do 
not apply any Galaxy plane cut-off and masks, nor do 
we dissect the whole sky into disjoint regions. Our work is 
closely related with Chiang & Coles (2000); Chiang (2001); 
Naselsky, Novikov & Silk (2002); Chiang, Coles & Nasel- 
sky (2002) and Chiang, Naselsky, Verkhodanov & Way 
(2003), from which some of the aspects of the CMB phase 
analysis were developed. A new element is to implement 
the phases of the maps from multifrequency K-W bands 
for reconstruction of the PCM. 

The layout of this paper is as follows. In Section 2 we 
describe the phases of the WMAP and illustrate the cross 
correlations between the foreground maps provided in the 
WMAP website and between the 5 frequency maps. We 
introduce the "non-blind" PCM method in Section 3 and 
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in Section 4 we introduce the "blind" PCM method and 
elaborate the 4 steps for foreground cleaning. We use a 
set of simulated maps, which we call the WMAP simulator 
as the numerical test on the "blind" PCM method. We 
produce the PCM using the "blind" method in Section 7 
and compare the PCM map and the power spectrum with 
those of the WMAP ILC map and the TOH foreground- 
cleaned map and Wiener-filtered map. We perform the 
Gaussianity test on these maps in Section 8 and the Con- 
clusion and Discussions are in Section 9. 

2. THE PHASES OF THE CMB AND FOREGROUNDS 

The statistical characterization of the temperature fluc- 
tuations of CMB plus foregrounds radiation and the in- 
strumental noise on a sphere can be expressed as a sum 
over spherical harmonics: 



AT(8^) = J2 a tm Y tm (0,tp), 



(1) 



t=0 m=- 



where ae m are the coefficients of the expansion. Homoge- 
neous and isotropic CMB Gaussian random fields (GRFs), 
as a result of the simplest inflation paradigm, possess 
Fourier modes whose real and imaginary parts are Gaus- 
sian and mutually independent. The statistical properties 
are then completely specified by its angular power spec- 
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In other words, from the Central Limit Theorem their 
phases 
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are randomly and uniformly distributed at the range 
[0, 2tt]. We hereafter denote at superscript "cm6" as the 
pure CMB signal. For the combined signal from the CMB 
and foregrounds we can write down 

|exp(i* £m ), (4) 
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where \ae m \ is the modulus and is the phase of each 
ai m harmonic. 

The statistical properties of the foreground signals are 
different from those of GRFs. In the contemporary CMB 
experiments the measured signals have contaminations 
from our Galaxy as well as from other extragalactic sources 
at all frequency bands. For example, for the K-Ka bands 
of the WMAP the foreground emissivity is mainly deter- 
mined by the combination of the synchrotron and free- 
free emissions. For higher frequency bands (Q, V and W 
bands) the contaminations include the free-free and dust 
components and some of the residues from the synchrotron 
radiation. Using different assumption on the properties of 
the foregrounds such as the frequency dependence for the 
amplitudes of the each kind of the foregrounds and detec- 
tor noise, the predicted spectral index for the point-like 
sources, statistical independence of the components ex- 
clude the Galaxy contaminations or the so-called "blind" 
Bayesian approach, several methods have been proposed 
for the separation of the pure CMB signal from various 
foreground contaminations (Tegmark & Efstathiou 1996; 
Bouchet & Gispert 1999; Hobson et al. 1998; Stolyarov et 

5 http://lambda.gsfc.nasa.gov 



al. 2002; Delabrouille, Patanchon & Audit 2003; Tegmark, 
de Oliveira-Costa & Hamilton 2003; Patanchon et al. 2003; 
Barreiro et al. 2003; Martinez- Gonzalez et al. 2003). 

Our PCM method, complementary to the methods men- 
tioned above, is based on a prevailing assumption: the 
phases of the "true" CMB signal should not correlate with 
those of the foregrounds. The advantage of this assump- 
tion is that it allows us without any special cut-offs of the 
Galactic plane to remove the foregrounds, in particular 
the Galaxy contaminations from the WMAP K-W maps. 
Below we use cross correlation of the phases between the 
WMAP bands to demonstrate pronounced cross correla- 
tions of the foreground phases for different bands, as shown 
in Fig 1. All the phases of the foreground components are 
derived from the WMAP foreground maps 5 . 
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Fig. 1. — The cross correlations of phases between all foreground 
maps for the K-W bands. The phases arc from all harmonic modes 
I < 50. All the foreground maps are taken from the WMAP website. 
Note also that on the bottom right panel we show the phases of the 
WMAP ILC map and the foreground map at the Ka band, which 
displays no correlations. 



In Fig. 2 we show the phase difference between the same 
?, m modes but for different K-W foreground maps, e.g. 



differences for Ka-Q band and Q-V bands are minimal, 
unlike, for example, for Ka-W bands. For comparison, in 
Fig. 3 we plot the cross correlations of phases for the total 
phases between the bands. As one can see from Fig. 1 and 
Fig. 3, all the properties of the foregrounds phases correla- 
tions manifest themselves in the cross correlation diagrams 
of the phases. 
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Fig. 2. — The phase difference at the same £, m modes for different 
K-W foregrounds maps, e.g. — against for the 

top left panel. 

In practice, we use the Healpix package (Gorski, Hivon 
& Wandelt 1999) to decompose each of the WMAP K- 
W maps for the coefficients of spherical harmonics 
where a = ai m and hereafter the superscript index j = 
1, 2, 3, 4, 5 correspond to the K, Ka, Q, V, W bands, respec- 
tively, after which we use the Glesp pixelization scheme 
(Doroshkcvich ct al. 2003) for the rest of data processing 
in order to control the accuracy of the harmonics. We 
obtain the phases = from all a^ (sec Eq. (2)- 
(4)). Each set of the a^ coefficients has the combination 
of different kinds of the foregrounds a^P'* (convolved with 
the beam) , the CMB signal & cmb (also convolved with the 
beam), and the instrumental noise n^: 
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where Efc 3 !"''' ^ + n^'- 1 = G^' is the common signal ex- 
cluding the CMB at the j-th band. Unlike the linearity 
in Eq.(5), the phases of the j-th band are related to the 
phases of each kind of the foregrounds, the CMB and the 
instrumental noise in a non-linear manner: 
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where £ denotes the phases of the CMB signal and 
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where H>g denotes the phases of the g-th foreground 
component at the j-th band and the phases of the 
noise. The phases of the common signal excluding the 
CMB for the j-th. band can be written as 
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Fig. 3. — The cross correlations of the phases between the maps 
of the K-W bands. 

As one can see from Eq.(7) and (8), the first term on 
the right hand side of Eq.(7) corresponds to the interfer- 
ence at the same multipoles £, m between different fore- 
ground components, whereas the second and the third 
correspond to between the noise-noise and between the 
foreground component-noise interferences at the same fre- 
quency band, respectively. We would like to point out that 
these interferences can be described in terms of the auto- 
and cross-correlations at the same channel. However, be- 
cause there is only one single realization for each frequency 
band these correlations are not well-established. It would 
be difficult to say, for example, that the instrumental noise 
does not correlate with the foregrounds tails by using only 
the 5 channels and the corresponding 5 values for each £, m 
harmonic. 
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In our phase analysis, therefore, we use some prelim- 
inary information about the possible values of the fore- 
ground and noise amplitudes at different multipole range. 
That is, we discuss the properties of the signals only at 
the range I < 50, so that we can avoid the beam shape 
and noise properties for all the K-W bands. This multi- 
pole range is also the available range for the preliminary 
foreground maps on the WMAP website. Moreover, we 
can assume that the contribution from the instrumental 
noise at this range is relatively small in comparison with 
the synchrotron, free-free and dust emissions (except for 
the K band). Generally speaking, we are investigating the 
properties of the CMB signal at the WMAP K-W bands 
at the multipole range comparable with the COBE range 
but slightly higher. At this range the contamination of the 
foregrounds is the major part of the signal, including the 
Galactic plane. 

3. THE "NON-BLIND" PHASE-CLEANED MAP (PCM) 

Note that we can neglect the instrumental noise and the 
beam shape properties when we consider multipole range 
I < 50. It is obvious that if the foreground signals are 
known, the pre-CMB signal can be obtained from a sim- 
ple subtraction of each band map by the corresponding 
foreground map. We call such derived map the pre-CMB 
signal because we assume that the foreground signals do 
not have the required accuracy. As a result some residues 
from the foregrounds propagate to the CMB signal. In 
reality, all properties of the foregrounds and the pre-CMB 
signal come from the same combination of the K-W maps 
and an uncertainties of the foreground and the CMB sig- 
nal estimation depend on each other. From a practical 
point of view, therefore, any "blind" methods would be 
applicable in order to find the correct phases of the CMB 
signal without any additional assumption about the statis- 
tical properties of the foregrounds (Tegmark & Efstathiou 
1996; Tegmark, de Oliveira-Costa & Hamilton 2003, for 
example). To make the first step for finding a "blind" 
method for the CMB reconstruction from the phases, we 
consider first a "non-blind" method and will suggest a sim- 
ple phase filter, the modification of which allows us to re- 
construct the pre-CMB phases by a combination of the 
ae, m coefficients of the initial CMB + foregrounds maps. 
This method generalizes the method from Tegmark & Ef- 
stathiou (1996), and the TOH approach, but in our case 
we require the minimal variance in the phases rather than 
in the derived map. 

We use all pairs of the WMAP maps for the reconstruc- 
tion of the pre-CMB maps. Below we use the pair of maps 
j = 1,2 (corresponding to the maps of the WMAP K 
and Ka bands) as an example to derive the weighting 
coefficients, it can nevertheless be applied to pairs such 
as j = 2,3 (Ka,Q), 2,4 (Ka,V) and 3,4 (Q,V). We ex- 
pect to obtain the derived map M from the pairs of the 
WMAP bands M« and M < 2 ) maps, whose spherical har- 
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Let us denote = |w^^ | exp(— itt^) where |w^| and 
are the moduli and the phases of the weighting co- 
efficients for the j-th band, respectively. Because of the 
normalization l) J+1 w^^ = I the moduli and phases 

of the filter are not independent 
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and the resultant phases of the derived map M are 
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where is the common signals excluding the CMB at 
the j-th band and is the phases of the j-th band. 

Thus, if we prescribe that the reconstructed phase \I> M 
from the weighting coefficients is that of the pre-CMB sig- 
nal (i.e. = £), we reach an exact solution for the 
moduli and phases for the weighting coefficients 
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|G« -G( 2 )| 
sinJ7 (1) = |w (2) |sinA, 
sinJ7 (2) = |w (1) |sinA, 



(13) 



where A = tl^ - fiW = *(2) _ 

Using the provided foreground maps from the 
WMAP database and their 5 maps for the whole signal 
we are able to reconstruct the pre-CMB signal using the 
filter Eq.(13). In Fig. 4 we plot the reconstructed images 
of the pre-CMB signal which contain some residues. For 
drawing all maps we use i < 50 only. The top panel is a 
simple subtraction of the V band map by the foreground 
map in the V band. The middle map was obtained us- 
ing the filter Eq.(13) and the pair of the WMAP Q and 
V maps and the their foreground maps for determination 
G (i) and A in the filter Eq.(13). The bottom panel is the 
difference between the first two panels. 

As one can see from Fig. 4 the derived pre-CMB maps 
look very similar, while the contaminations from the bright 
point-like source and diffuse residues exist notably along 
and perpendicular to the Galactic plane, respectively. In 
order to compare our derived PCM and the WMAP maps 
we plot the power spectrum of the 3 maps in Fig. 5, using 
a simple pseudo-CV estimator (Hivon et al. 2002), 
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Note that our PCM has smaller power for the whole 
range of the multipoles in comparison with that derived 
from the V band map. The flatness of the power spec- 
trum at I > 20 clearly demonstrates the contamination of 
the bright point-like source residues mainly concentrated 
along the Galactic plane. 
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Fig. 4. — The maps for the derived pre-CMB signals. Top panel 
is the map taken from a direct subtraction of the foreground in V 
band by the WMAP total signal in V band. The middle panel is the 
cleaned-up map obtained from the "non-blind" PCM method from 
the WMAP Q and V bands, and the bottom the difference between 
the first two panels. 
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4. THE "BLIND" RECONSTRUCTION OF THE CMB 

In this section we describe the "blind" PCM method 
for separation of the CMB signal from the foregrounds us- 
ing the K-W band signals at the harmonic range £ < 50, 
avoiding again the beam convolution and the instrumental 
noise in Eq.(5)-(8). The main assumption is that in the 
different frequency bands the phases of the foregrounds do 
not correlate with those of the true CMB signal. The re- 
sult of separation improves if there are strong correlations 
between the phases of the maps between the compared 
frequency bands (see Fig. 3). 

The basic idea of the "blind" PCM method is to gen- 
eralize the minimization scheme by TOH and Tegmark & 
Efstathiou (1996), including the minimization of the cross- 
correlations between derived CMB signal and foregrounds 
phases. This "blind" method consists the following four 
steps and does not need any galactic cut-offs and dissec- 
tion of the whole sky into disjoint regions. 

a. We firstly group the following combinations 
of the WMAP maps: Ka-Q, Ka V and Q V. They 
are grouped due to highly correlated phases between the 
bands shown in Fig. 3. We did not include the pair of the 
W and K bands in the separation procedure because of 
the reasons mentioned in Appendix B. 

b. For each pair of the maps (Ka-Q, Ka — V 
and Q — V) we want to reach the pre-CMB map 
a^L = cim derived from the pair of the bands. For 

Ka and Q (j — 2, 3), for example, we set = a£2 an( l 

M^> = af^ an d use optimization coefficients = 
similar to the TOH and Tegmark & Efstathiou (1996): 

3=2 3=2 

The coefficients are subject to the following con- 

straints 
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where S/ST^' are functional derivatives. The then 
have the following forms 
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where are the phases of the total signal at map. 6 
After the pre-CMB separation, we obtain the pre-CMB 
signal C£ m and the pre- foreground at each bands. 

c. The next step is to find the filter g^J to im- 
prove the foreground separation using the phase 
variance minimization. We will use the filter in 
Eq.(19), which minimizes the phase difference between re- 
constructed and the CMB phases. We assume that g^' 

6 We want to emphasize that one may think it could be possible to omit the summation other m in Eq. (16) and Eq. (17) in order to optimize 
the phases at each I, m mode. Then r' J ) depend on m too. The reasons why we use the filter with the form Eq.(17) are given in Appendix A. 



Fig. 5. — The power spectrum for the derived from the WMAP V 
band map (the dash dotted line), the PCM (the solid line) and the 
difference between the two maps (the dash line). 
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does not possess any imaginary part and depends only on 
I but not on m. The theoretical basis of our choice of the 
form of the filter is for minimization of the weighting phase 
variance 

i \ n cmb\2 

V=— -V l -^L(*-n 2 ^min (18) 

This condition is equivalent to minimization of the error 
e = {21+ l)- 1 J2 m \af m - a c t,™ b \ 2 between reconstructed 
{af m ) and the true (a^ b ) CMB signals 7 . Using fore- 
grounds estimator G l = a* — a M from the step b, we 
obtain the following filter gy> instead of the filter (13). 

(2) _ T. m |G< 3 >|[|G( 3 »|-|G( 2 »|cos(*( 3 »'-f-*< 2 »-')] 

5 — £ m G( 2 >-G< 3 >| 2 

(3) _ S: m |G' 2) l[|G( 2 )|-|G( 3 )|co S (^).^^)^)] 

6 — -£ m G( 2 >-G< 3 >| 2 ^ ^> 

where \|>(.j)>/ are the phases of the foregrounds at the j-th 
band. This filter provides a new estimation of the fore- 
grounds and the pre-CMB signals which we shall use for 
the next step of the subsequent iterations by the filter 
Eq.(19). So, at each iteration we insert the foreground 
moduli and phases obtained from the previous iteration 
to the same form of the filter Eq.(19). In general, 3 such 
iterations are enough to produce a convergent map. 

d. We use a MIN-MAX filter to further clean up 
the residues. After the step c. we have three cleaned- 
up maps from the WMAP map pairs: Q-V, Ka-V and 
Ka-Q bands. These 3 maps correspond to the minima of 
the phase difference (see Fig. 3). They nevertheless still 
contain some of the foreground residues, notably along 
the Galactic plane. In order to detect and minimize such 
a residues we used a MIN-MAX filter on the same pixel 
among the 3 cleaned-up maps. The definition of the MIN- 
MAX filter is the following. Using maps with the same 
pixelization scheme, we can designate all the pixels in the 
same way and compare simultaneously at the same pixel 

p these 3 cleaned-up maps ATp^ by pairing them up, 
where the index corresponds to the 3 cleaned-up maps. 
The MIN filter L m i n allows us to find the signal which 
has a minimal (absolute) amplitude at the pixel p, i.e. 

L min ({AT^}) -» AT™ in , where |AT™ in | = min{|AT p (i) |}. 

In other words, we take the smallest absolute amplitude 
at the pixel p from the cleaned-up maps and attach its sign 
for that pixel (or even more concisely, the value which is 
the closest to zero) . 

The reason for such a filter is clear. The signal at each 
pixel is a combination of the underlying CMB signal and a 
small residue from the foreground cleaning. If the derived 
CMB signal correlates little with the foregrounds, it is nat- 
ural to expect that all deviations of the ATp^ at p of the 
different maps are caused by the residues but not by the 
CMB. In addition to the L m i n filter we can define a MAX 
filter, L max ({AT p W }) ->■ AT™ ax which is the farthest value 
from zero between the maps at pixel p. This MAX filter 
reproduces the most possible contribution from the fore- 
ground residues at each pixel. Thus, we find our final 
PCM, the true CMB map according to the criterion that 
AT p maK — AT™ m map has the mininum power. We apply 

7 See Appendix B for details 



our MIN-MAX filter on Q-V, Ka-V and Ka-Q cleaned-up 
maps and obtain the PCM signal from Q-V, Ka-V maps 
combination. 

At the end of this Section we note the following com- 
ment. In Appendix C we show analytically that if the 
phase £ of some ae m of the CMB signal is close to ir/2 
or 37r/2 the error of the phase reconstruction is very high. 
Practically these phases can not be reconstructed. We call 
this problem the ll 7r/2" problem. Our numerical tests (see 
section 5) confirm this conclusion. 




Fig. 6. — The reconstructed PCM from the WMAP simulator. Top 
panel is the simulated CMB map and the middle is the reconstructed 
PCM from the two pairs of simulated bands, Q— V and Ka-V. The 
bottom panel is the difference between the two. 



5. THE NUMERICAL TEST 

To numerically test the "blind" PCM method we use 
a WMAP simulator to mimic all the properties of the 
WMAP foregrounds at K-W bands. We again exclude the 
instrumental noise and the beam shape properties, which 
are insignificant for the multipole range I < 50. We take 
the best-fit WMAP ACDM power spectrum to produce a 
simulated CMB sky with random phases. Such CMB sky 
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is truly a Gaussian random realization. The morphology 
of this simulated CMB map is obviously different from 
that of the WMAP ILC and the TOH maps (since it has a 
different set of random phases). For our wmap simulator 
we adopt all the foreground maps for each frequency band 
from the WMAP website. We also confirm that the simu- 
lated CMB signal does not have any cross-correlation with 
the foregrounds and it has well-defined statistical proper- 
ties (the power spectrum and random phases). 
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Fig. 7. — The reconstructed power spectrum of the WMAP sim- 
ulator from the "blind" PCM method. The solid line is the power 
spectrum of the simulated CMB signal and the dash-dot line is that 
of the reconstructed PCM. Note that both curves almost overlap 
except for £ = 2,5 and 7. The dashed line is the difference of these 
two. These 3 power spectra correspond to the 3 panels in Fig. 6. 
The reconstructed PCM signal is produced from two pairs of the 
simulated channels, Q-V and Ka-V. 
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Fig. 8. — The cross correlation of the phases between the sim- 
ulated (random-phase) CMB and the reconstructed CMB (PCM) 
signal. The solid diagonal line indicates an asymptotic when the 
simulated and reconstructed phases are identical. 

In Fig. 6 we show the result of the CMB reconstruction 
after 10 iterations from the "blind" PCM method. As 
we mention above the third and all subsequent iterations 
practically produce a convergent cleaned-up map. We take 
two pairs of simulated channels, Q-V and Ka-V, to recon- 
struct the signal. The result of reconstruction for both 
pairs was practically the same, and the MIN-MAX filter is 
used to produce the resultant reconstructed map shown in 
Fig. 6. In Fig. 7 we plot the power spectra for the simulated 
and the reconstructed CMB (PCM) signals as a function 
of I. Both curves almost overlap except for t = 2, 5 and 
7. 

In Fig. 8 we plot the cross-correlation of the phases be- 
tween the simulated CMB signal and the reconstructed 
CMB (PCM). As one can see from Fig. 8 for reconstructed 



CMB signal the quadrupole component has smaller power 
than the simulated CMB. The reason for that is that the 
peculiar phases of the quadrupole component Q 2 o, Q21 
and Q 2 2- Accidentally, in our simulated realization of the 
CMB (see Fig. 6) the phases of the component Q20 and 
Q21 are zero whereas for the Q22 component the corre- 
sponding phase is close to 3n/2 with only 10% deviation. 
That is why the reconstructed phase of Q22 component 
has £22 = 6.227 radians. 

In Fig. 9 we plot the cross correlation of the phases be- 
tween the reconstructed CMB signal (PCM) and the sim- 
ulator foregrounds at all 5 simulated bands. The scat- 
tering of points in all 5 panels indicates that the phases 
of the PCM display no serious correlations with those of 
the foregrounds. The WMAP simulator therefore confirms 
the usefulness of our "blind" PCM method for the CMB 
extraction without any assumptions about the statistical 
properties of the foregrounds. 
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Fig. 9. — The cross correlations of the phases between the recon- 
structed CMB signal (PCM) and the simulator foregrounds of all 5 
simulated bands. The scattering of points in all 5 panels indicates 
that the phases of the PCM display no correlations with those of 
the foregrounds. 



6. THE PCM DERIVED FROM THE 1-YEAR WMAP DATA 

For the reconstruction of the CMB signal at I < 50 
from the K-W bands of the WMAP data we use the three 
pairs of the maps: Ka-Q, Q-V and Ka-V. We repeat 
the steps b and c of the "blind" PCM method : each 
of them has been convolved by the weighting coefficients, 
Eq.(17) and then by the coefficients Eq.(19) using only 2 
iterations. All the next iterations do not change signifi- 
cantly the a,f m coefficients (the corresponding error being 
less than 10~ 3 %). The results are present in Fig. 10, in 
which all the 3 cleaned-up maps have pronounced simi- 
larity in morphology, albeit with some residues from the 
galactic foregrounds. 
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Fig. 10. — The cleaned-up maps after the steps b and c from 
the "blind" PCM method on the 3 pairs of WMAP maps. The top 
panel is the cleaned-up map from the Ka-Q maps, the middle the 
cleaned-up map from the Ka-V maps, and the bottom from the 
Q— V maps. 

In Fig. 11 we show the result after the step d, the MIN- 
MAX filtering. The top panel is our CMB signal (the 
PCM), the middle is the WMAP ILC map for comparison 
and the bottom is the difference of the two. 

In Fig. 12 we plot the CMB angular power spectra from 
different methods: the thick solid line is the PCM, the 
dashed line the WMAP ILC map, the dash-dot line the 
TOH foreground-cleaned map, and the thin solid line the 
TOH Wiener-filtered map. The PCM power spectrum is 
smaller then ILC power spectrum at the whole multipole 
range £ < 50 and it reproduces a similar power spectrum 
as the TOH Winer-filtered map. 

In Fig. 13 we show the cross correlation of the phases 
between the PCM and the foreground maps of the 
WMAP bands. Again the scattering of the points shows 
that the phases of the PCM practically have no correla- 
tions with those of the foregrounds. 




Fig. 11.— The CMB map (the PCM) after the step d of the 
"blind" PCM method (top panel). The MIN-MAX filter is per- 
formed on the 3 panels of Fig. 10. For comparison, the middle is the 
WMAP ILC map and the bottom is the difference of the two. 
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Fig. 12. — The comparison of the CMB angular power spectra 
from different methods. The power spectrum of the PCM is shown 
with the thick solid line, that of the WMAP ILC map is shown with 
the dashed line, the TOH foreground-cleaned map the dash-dot line, 
and the TOH Wiener-filtered map the thin solid line. 




Fig. 13. — Cross correlation of the phases between the PCM and 
the WMAP foreground maps of the WMAP Ka, Q, V and W bands. 
Again the scattering of the points shows that the phases of the PCM 
practically have no correlations with those of the foregrounds. 
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Fig. 14. — Cross correlation of the phases between the PCM and 
the WMAP ILC map (top left), between the PCM and the TOH 
foreground-cleaned map (top right), and between the PCM and the 
TOH Wiener-filtered map (bottom). 

Although the phases of the PCM have strong correla- 
tions with those of the ILC map, the TOH foreground- 
cleaned map and the Wiener-filtered map, as shown in 
Fig. 14, there are undoubtedly deviations in phases be- 
tween them, which is an indication of different morphology. 

The CMB map from the "blind" PCM method allows 
us to reconstruct the foregrounds for each WMAP band. 
In Fig. 15 we show at the top the foreground map at 
the V band, a simple subtraction of the PCM from the 
WMAP V-band map. In the middle panel we show the 
WMAP foreground map at the same band and the dif- 
ference of these two panels is shown at the bottom. 
As one can see from this Figure, an additional power 
of the foregrounds from the PCM subtraction than the 
WMAP foreground map is clearly shown in the third 
panel, in particular we can see strong contribution of the 
of the point-like sources residues along the Galactic plane 
which are detected and removed as foregrounds from the 
PCM subtraction. 
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Fig. 15. — The foreground map at the V band from the sub- 
traction of the PCM from the WMAP V-band map (top), the 
WMAP foreground map at the V band (middle) and the difference 
of the first two maps (bottom). 



7. NON-GAUSSIANITY TESTING ON THE PCM 

In this section we use the phase mapping technique 
as a non-Gaussian test to compare the non-Gaussianity 
between the PCM, the WMAP ILC map and the TOH 
Wiener-filtered map. This method has been successful in 
detecting non-Gaussianity in the derived TOH maps from 
the WMAP data (Chiang, Naselsky, Verkhodanov & Way 
2003). 

The null hypothesis of the phase mapping technique is 
based on the so-called random-phase hypothesis, which 
is served as a practical definition of homogeneous and 
isotropic Gaussian random fields (Bardeen et al. 1986). 
Therefore, any non-random phases from harmonic anal- 
yses is interpreted as manifestation of non-Gaussianity. 
Due to the circular character of phases, we use the return 
mapping of phases with certain separations (Am, A£) to 
render phases on return maps. Testing random phases 
is thus equivalent to testing randomness of each return 
map. There are many ways of quantifying such ensemble 
of return maps. We choose the simplest x 2 statistic as in 
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Chiang, Naselsky & Coles (2002): we first use a Gaussian 
window function to smooth each return map and use a sim- 
ple x 2 = (V-W) j [p(*j j) ~ p] 2 /p to yield the statistics, 
where p is the mean value and M is the number of pixels. 
The null hypothesis will have mean (x 2 )p = l/47ri? 2 and 
dispersion S| = l/n 3 R 2 (M/2) . 

We perform the mean x 2 statistics derived from phase 
mapping technique for the PCM, the ILC map and the 
TOH Wiener-filtered map. We take the phases from £ = 2 
to maximum £ = 40 and render 256 points on each re- 
turn map. The smoothing scale of the Gaussian window 
function is R = 1. The ensemble of return maps consists 
of separation Am from —9 to 20 and A£ from 1 to 15. 
From the distribution function of such ensemble, we use 
the mean (x 2 ) and the dispersion S of the distribution 
function as two independent variables for displaying the 
statistics. We simulate 2000 realizations of Gaussian ran- 
dom fields to yields reliability. In Fig. 16 the dashed and 
dotted contours indicate 68%, 95% confidence level for de- 
tection, respectively. The non-Gaussianity through phase 
mapping test shows that the PCM and the Wiener-filtered 
map have close non-Gaussianity, whereas the WMAP ILC 
map shows less non-Gaussian. Note that, however, both 
the ILC map and Wiener-filtered map are produced by 
dissecting the whole-sky CMB maps into 12 and 9 dis- 
joint regions for cleaning, which systematically change the 
phase configuration for the whole sky. 
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Fig. 16. — Non-Gaussianity test using the phase mapping tech- 
nique for the PCM map (cross) at multipolc range 2 < I < 40, 
the WMAP ILC map (triangle) and the TOH Wiener-filtered map 
(box). The dashed contour corresponds to 68% CL, the dotted con- 
tour corresponds to 95% CL. 

8. DISCUSSIONS AND CONCLUSIONS 

In this paper we have proposed a new method for 
the separation of the CMB signal from the Galactic and 
other foregrounds based on the phase cleaning of the 
WMAP K-W bands at the multipole range £ < 50. This 
method is based on a simple but prevailing assumption 
that the phases of the CMB signal should not correlate 
with those of the foregrounds. We first propose a "non- 
blind" PCM method using the foreground maps available 
at the WMAP website as preliminary foregrounds. We 
then introduce a "blind" PCM method with 4-step fore- 



ground cleaning. We firstly group into pairs from the maps 
which have strong cross correlations in phase. We then in 
step b use a simple TOH minimization of the variance of 
the combination between different pairs of K-W maps. In 
the step c we propose a phase-cleaning filter for each de- 
rived map from the constraint of minimal variance of phase 
difference. Finally, we propose a MIN-MAX filter in order 
to obtain the PCM map. Our CMB signal (the PCM) has 
the power spectrum in an excellent agreement with the 
best fit WMAP cosmological model, but the phases are 
slightly different from those of the WMAP ILC map. We 
also confirm that, after phase cleaning, the PCM has essen- 
tially no correlation in phases with their cross-correlation 
with the foregrounds. We would like to point out that the 
PCM method can be used with different component sep- 
aration methods such as the Maximum Entropy Method 
(MEM) or another methods, which can be used an initial 
step before the phase cleaning. Such generalization on the 
PCM method is in progress. 

Using the PCM method we have found corresponding 
corrections of the foreground maps for each K-W band, 
which is useful for the data analysis for the upcoming 
Planck mission. The power spectra of the common fore- 
ground maps from our PCM subtraction arc slightly dif- 
ferent from the WMAP foreground maps. Moreover, the 
PCM does not show any serious contamination from the 
point-like source residues mainly thanks to the MIN-MAX 
filter. 

The PCM CMB signal, as it follows from comparison 
between our CMB and WMAP maps has a power smaller 
then ILC map. 

We also confirm that the PCM displays a low level 
of power at the quadrupolc component, which also 
manifest itself in the WMAP ILC map and the TOH 
map (Tegmark, do Oliveira-Costa & Hamilton 2003; de 
Oliveira-Costa, Tegmark, Zaldarriaga & Hamilton 2003). 
However, the low level in C'e is typical at the multipolc 
range 2 < £ < 10. In order to compare the power spec- 
trum of the PCM with that of the WMAP best fit model 
and the TOH Wiener-filtered map, in Fig. 17 we plot these 
power spectra with the error bars caused by the cosmic 
variance. 
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Fig. 17. — The power spectra of the PCM (the thick solid line) 
and the TOH Wiener-filtered map (thick dashed line). The cosmic 
variance error bars with 68% CL arc indicated with crosses with the 
triple dots dashed lines. 

As one can see from Fig. 17 all deviations between the 
PCM and the TOH Wiener-filtered map are within the 
68% CL. Some of the multipolc range, however, seems 
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to be quite peculiar. Namely, at 20 < £ < 24 the power 
of both the PCM and the TOH Wiener-filtered map has 
some suppression. Does it indicate that we have a pe- 
culiar realization of a Gaussian random CMB field? Or 
we still have the contaminations of the foregrounds which 
mimic as extra CMB power? All these issues can only be 
discussed when more data are available, particularly the 
future polarization observational data. 

In this paper we concentrate only on the low multipole 
range of the CMB £ < 50, at which the beam shape and 
the instrumental noise are not crucial. The development 
of the PCM method and the results at the multipole range 



£ > 50 is in preparation. 
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APPENDIX 



A : THE CONCEPT OF THE "BLIND PCM METHOD 

From Eq.(9) and (13), one might think that it is possible to generalize the "non-blind" phase reconstruction method 
to the "blind" variant simply using the ag m coefficients in Eq.(13) as an estimator for the reconstructed phases of the 
pre-CMB signal. Such generalization, however, would be incorrect because of the following reason. Below we consider 
only for the multipole range I < 50 to avoid the instrumental noise and beam shape properties. According to Eq.(5) the 
moduli of the ae m coefficients are related with the moduli of the foregrounds and CMB for each WMAP band as (with 
the phases in Eq.(5) being replaced by tyW'f , the resultant foreground phases) 



|G (i) | 2 + 2|G^)||a cmfc |cos(* 0%/ 
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« |G^|. Using Taylor series one can obtain 
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£). Then, from Eq.(12) we 
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(A2) 



The comparison of Eq.(A2) and Eq.(12) allow us to conclude that the deviation of the reconstructed phases from the 
CMB phases should be very high even for the small values of A-parameter. This means that the reconstructed phases 
of the CMB signal would be close to the foreground phases and the corresponding power spectrum should have very 
significant errors. It is important to note that such a conclusion means that in order to extract the correct phases 
and amplitudes of the CMB signal from the data using a "blind" method of separation, it is necessary to decrease the 
contamination in Eq.(Al) of the linear term (~ |a cmh | cosfifW'J — £)). The only way to do this is to take into account 
that J^m |a cmb |c°s(*°' ) ' / - -> for hi g h e > m - In such a case £j(-l) j+1 wW|a W l ~ C(|a cmb | 2 ) and corresponding 

8 http : //www . eso . org/science/healpix/ 
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error of the CMB phase reconstruction should be of the order la™ 16 ! / min {|G^|}. Moreover, averaging over the m values 
means that the corresponding weighting coefficients is a function of I only, but not of m, which is ideologically close 
to Tegmark & Efstathiou (1996) and the TOH method, where = w^\i). We would like to point out, however, that 
for such an optimization, as indicated in Eq.(12), the reconstructed CMB phases will have the cross correlations with the 
foregrounds phases at different levels for different modes. 



B : MINIMIZATION OF THE PHASE DIFFERENCE 

Assuming that the reconstructed phase should be close to the CMB phase £, from Eq.(12), we obtain 

V. gO") |GW) I sin(*«^ - cos * M 



arcsin 



£\gO')|GC?')|cos*OV + |a cmb |cos£' 



(Bl) 



As mentioned above, a new element of the "blind" reconstruction of the CMB signal is the correlation between the derived 

linear combinations 

g(j)|GW)|cos*^, (B2) 



phases \I/ M and the foreground phases Vt^'^, ty^'* f° r au linear combinations of the maps. Moreover, if 



|a cmfc |cos£»E' 

3 



then \I> M ~ £ while for some of the £,m modes which corresponds to inequality |a cmh |cos£ <C ■ | cos , 

reconstruction of the CMB phases needs more specific investigation, which will be discussed at the end of the section. 
For all £, m modes satisfying Eq.(B2) in order to decrease the correlation between the foregrounds and CMB phases we 
introduce the optimization scheme based on the iterative solution of the Eq.(Bl). The first step of iterations is to minimize 
the variance 

/.2-7T 
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which corresponds to the equation 
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Note that /? could be both positive and negative depending on the phases i&W'f . As we mention above, Eq.(B4) has an 
exact solution not for all values of £, m. For £, m satisfying Eq.(B2) we can neglect the (3 parameter in Eq.(B5) and obtain 
h ~ 2tt, I 2 ~ 0. Thus, from Eq.(B4)-(B6) we have 



g (D = 
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(B7) 



The filter Eq.(B7) depends on both I and m. As one can see from Eq.(B7), if the phases of the foregrounds \I>( 2 )'.f 
and vl^ 1 )*/ are identical, then one can have the same form for the g^ 1 ) and g( 2 ) weighting coefficients as in Eq.(13). 
Unfortunately, because of the non-zero phase difference \J>( 2 )^ — ifrWJ the filter (B7) is not useful for the "blind" method 
of the CMB reconstruction, ft is due to the instability of the iteration sheme, which reproduces at the very first step the 
phases of the foregrounds and use them as an attractor. 

Now following the concept in Appendix A, let us introduce another concept of the phase optimization based on the 
minimization of the weighting phase variance 
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which corresponds to the following 
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It is easy to demonstrate that the filter Eq.(BlO) corresponds to the minimum of the variance ^ m |a M — a cmb \ 2 . We 
propose to use this filter for optimization of the pre-CMB reconstruction at the second and subsequent iterations. 



C : THE "7r/2" PROBLEM OF THE PHASES RECONSTRUCTION 

In this Appendix we discuss some problems of the CMB image reconstruction related with optimization methods. The 
exact solution for the reconstructed phases in a form of Eq.(BlO) has some peculiar reconstructed phases \I> M if the CMB 
phases £ are close to tt/2 and 37r/2. In such a case, from Eq.(B4) we obtain 

^ ( \a- b \cosS E jg (j) |G (j) |sin^0-)-/ \ 

where £ = 7r/2 + S, S -C 1 and |a cm(> |cos£ -C g^\G^\cos^/^^ . We point out that depending on the foreground 

phases the reconstructed phase ^ M could be arbitrary. It means that phases close to tt/2 will be reconstructed with a 
big error. The same argument can be applied to the case when £ is close to 37r/2. 



